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It is commonly assumed that the cosmological and astrophysical bounds on the mixings of sterile 
with active neutrinos are much more stringent than those obtained from laboratory measurements. 
We point out that in scenarios with a very low reheating temperature Trh << 100 MeV at the end 
of (the last episode of) inflation or entropy creation, the abundance of sterile neutrinos becomes 
largely suppressed with respect to that obtained within the standard framework. Thus, in this case 
cosmological bounds become much less stringent than usually assumed, allowing sterile neutrinos 
to be "visible" in future experiments. Here, we concentrate on massive (mostly sterile) neutrinos 
heavier than 1 MeV. 

PACS numbers: 14.60.St, 98.80.Cq 

I. INTRODUCTION 

In inflationary models, the beginning of the radiation dominated era of the Universe results from the decay of 
coherent oscillations of a scalar field and the subsequent thermalization of the decay products into a thermal bath 
with the so called "reheating temperature", Trh. The standard computation of the relic densities relies on the 
assumptions that Trh was large enough for the particles of interest to have reached thermal equilibrium and that 
the entropy of matter and radiation is conserved after they decouple. However, there are non-standard cosmological 
models in which these assumptions about the epoch of the Universe before Big Bang Nucleosynthesis (BBN), an 
J> ■ epoch from which we have no data, do not hold. These include models with moduli decay Q-ball decay @, and 
thermal inflation [3]. In all of these models there is a late episode of entropy production in which the Universe is 
reheated to a low Trh. This temperature may have been as low as ~4 MeV [4[ once all cosmological data are taken 
into account, while BBN-only data impose a lower bound of ~2 MeV on Trh [E[ if active neutrino oscillations are 
taken into account (up from ^0.7 MeV [f| if they are not). It is well known that a low reheating temperature inhibits 
the production of particles which are non-relativistic or decoupled at T & Trh @, H, @ ■ The final number density of 
active neutrinos starts departing from the standard number for Trh ^ 8 MeV but stays within 10%, 20% or 50% of it 
for Trh 5 MeV, 4 MeV and 3 MeV, respectively. For Trh = 1 MeV the number of tau- and muon-neutrinos would 
be about 2.7% of the standard number [a, 1 1 01] - 

Low Trh cosmological scenarios are more complicated than the standard ones. Different aspects of these models 
have been studied with interesting results, but no consistent all-encompassing scenario exists yet. Baryogenesis could 
be produced through the Affleck-Dine mechanism, in a model similar to that of Ref. Dark matter could still 

consist of the lightest supersymmetric particle or other Weakly Interactive Massive Particles (WIMPs), produced 
either thermally or non-thermally [l2j • Alternatively, MeV scalars have been proposed as DM candidates [H, : 
they would be in thermal equilibrium in the Early Universe with neutrinos, electrons, positrons and photons and 
would freeze-out at 100 's keV- MeV temperatures. 

If sterile neutrinos exist and have no extra standard-model interactions, the dominant mechanism of production in 
the early Universe is through their mixing with active neutrinos (l5| . Dodelson and Widrow [lj| (see also Refs. fl7lll8l]) 
provided the first analytical calculation of the production of sterile neutrinos lighter than 100 keV in the early Universe, 
under the assumption (which we maintain here) of a negligible primordial lepton asymmetry. 

In general, given that the production rate of sterile neutrinos of mass m s is maximum at a temperature T max ~ 
1.3 GcV (mg/MeV) 1 / 3 if the reheating temperature is smaller, Trh < T max , the production of sterile 

neutrinos is suppressed. Hence, the main idea of this paper is that the primordial abundance of sterile neutrinos 
does not necessarily impose their mixing with active neutrinos to be small. A low reheating temperature scenario 
would suppress the sterile neutrino production, weakening the cosmological bounds. Thus, it might be possible to 
consider massive (mostly sterile) neutrinos of any mass and mixing with active ones, as long as laboratory bounds are 
satisfied. These neutrinos could, therefore, be revealed in future experiments. We concentrate here on the production 
of massive (mostly sterile) neutrinos heavier than 1 MeV through the conversion of active neutrinos for Trh < m s , 



OO 

o 



- 1—1 

X 



2 



having already applied the same ideas to lighter sterile neutrinos 0]. By using different approximations, we obtain 
an analytical result for the sterile neutrino abundance. In this way, we are able to write all our results in a simplified 
form. Although it is a rough approximation, it allows us to have a qualitative understanding of the problem. For 
lighter neutrinos, m s < 1 MeV, the analytical calculation of Ref. [9j of the final sterile neutrino abundance turned 
out to be correct within an order of magnitude [icil ]. Like in Ref. [J, here the active neutrinos are assumed to have 
the usual thermal equilibrium distribution f A = (exp(E/T) + l)" 1 with E = p. Thus, following Ref. 0, we restrict 
ourselves to reheating temperatures Trh > 4 MeV. 

The paper is organized as follows. In Section |TT] we obtain an analytical approximation for the sterile neutrino 
abundance as a function of the reheating temperature. In Section IIIII we describe the laboratory bounds for sterile- 
active neutrino mixing for each of the three active flavors. Using the result of Section [III we compute in Section IIVI 
the astrophysical and cosmological bounds and show how, in some cases, they are completely evaded. Finally, in 
Section El we draw our conclusions. 



II. STERILE NEUTRINO ABUNDANCE 



For simplicity, our analysis is based on the two-neutrino mixing approximation. In this way, an analytical un- 
derstanding of the problem is possible. Within this approximation, the vacuum mixing angle sin# represents the 
amplitude of the heavy mass eigenstate v 2 in the composition of the active neutrino flavor eigenstate v ai i.e., 
v a = cos 8 v\ + sin# 1/2, v s = — sin# v\ + cosf? for a — e, fi, r, where v\ is the light mass eigenstate and v 2 
is the heavy mass eigenstate, which for small sinf? is mostly sterile and whose mass we call m s . 

In order to obtain the distribution function of sterile neutrinos at a given temperature after the last episode of 
inflation, we start from the Boltzmann equation , 

--Hp—)f s =I con , (2.1) 



dt dp 

where 7 co u ~ r s (T) (/° q — f s ) is the collision integral. Here, /° q = (exp (E/T) + 1) _1 is the Fermi-Dirac distribution 
that heavy neutrinos would have if they were in thermal equilibrium. In many cases the approximation T oc 1/a, 
with a the scale factor, is sufficiently accurate. Hence, by using T — —HT, Eq. (|2.1[) can be rewritten as (see e.g. 
Ref. [H), 

- HT (^k) ^r s (T)(/ s ^-/ s ). (2.2) 



dT J 



p/T 



For neutrino masses much smaller than the temperature of the plasma, (p) ~ (E), (m s < 1 MeV), the averaged 
rate of sterile neutrino interactions is given by 

r s (T) ~ ~ sin 2 26 nl d a G 2 F ET 4 , (2.3) 

where 9 m is the mixing angle in matter, Gf is the Fermi constant and d a — 1.13 for sterile neutrino mixing with 
v a = v e and d a = 0.79 with v a = i^, T . For T < 1.5 GeV (rris/MeV) 1 ^ matter effects are negligible and hence 
sin 2 26 m ~ sin 2 26 is a very good approximation. Plugging Eq. (|2.3[) into Eq. (|2.2I) . and solving for f s in the limit 
f s << /g q , the distribution function of massive (mostly sterile) neutrinos lighter than 1 MeV was found to be @ 

f s (E,T) * 3.2 d a ( ^) 3 sin 2 20 (|) , (2.4) 

for Trh *C T max . This distribution results in a number density of light (m s < 1 MeV) sterile neutrinos given by 

n s ~ 10 d a sin 2 29 ( n a , (2.5) 



where n a = 0.09 g T 3 is the number density of a relativistic fermion with g degrees of freedom in thermal equilibrium. 
Notice that the number density of sterile neutrinos depends on both the active-sterile mixing angle and the reheating 



1 The distribution functions are defined for mass eigenstates (see the Appendix). 
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temperature. A low reheating temperature implies a small sterile number density, even for active-sterile mixing angles 



as large as experimental bounds permit (see below). However, Eq. 



is only valid if the condition f s « /| q is 



satisfied. We would like to extend here this result, for the case when this approximation starts to fail. This was not 
included in Ref. 0]. Nevertheless, Eq. (|2.2p can be solved perturbatively and Eq. (|2.4p should be replaced by 



f s {E,T)~ (1 - e~ s ) f°*(E,T) 



(2.6) 



where S = 3.2 d a (Tr/5 MeV) 3 (E/T) sin 29 is the coefficient multiplying the equilibrium neutrino distribution in 
Eq. (|2.4[) . The sterile neutrino number density n s results from a numerical integration of this distribution. This 
perturbative solution is valid for S < 1. 

For heavier sterile neutrinos, with m s > 1 MeV, and for the range of temperatures explored here, the heavy neutrino 
mass needs to be taken into account and the averaged production rate of sterile neutrinos T s is given by (l9l . |20| 



r.CO = - 
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(2.7) 



where the first term is due to inverse decay and the other terms correspond to two-to-two particle processes. The last 
term in parenthesis, ~ T 4 /T s m 5 s , due to oscillations, is the only one remaining as m s —> 2 . The function r s in the 
denominator is the heavy neutrino lifetime. As in the case of m s < 1 MeV, Eq. (12. 7| is valid when matter effects are 
not important. 

For m s < 77V ~ 140 MeV, the massive (mostly sterile) neutrino can decay into a light neutrino and two leptons 
v s — > v a + I + I, mainly v a vv and v a e + e~ . If the active neutrino mixing with the sterile is v T or v^, the decay of v-2 
happens through neutral currents and the lifetime is 



1.0 sec /10 MeVV 



sin 2 26> 
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If instead v s mixes mostly with v e , the factor 1.0 sec should be replaced by 0.7 sec [19|, |20|, due to the presence of 
charged currents. However, we are not going to keep this distinction in the following. For < mi < 2rn M , j/2 
decays mostly into n + e~ and n~e + , and the decay is much faster than Eq. (|2.8|) JJ,[2l[. For even larger masses, 
other decay modes open up. In the following we will restrict ourselves to the range m s < m n ~ 140 MeV, which is 
enough for our purposes of showing the main characteristics of the low reheating temperature cosmological scenarios 
we envision. 

Assuming that the bulk of sterile neutrinos are produced after the reheating of the Universe 3 , namely that f s ~ 
at T — Trh, we solve analytically Eq. (|2.2p for Trh < m s , after plugging Eqs. (|2.7p and (|2.8p into it. In order to 
analytically solve the equation, we will make several approximations. First, we assume that the actual distribution 
function of the heavy neutrinos is always much smaller than the equilibrium distribution, f s << f° q . Then, we 
approximate the Fermi-Dirac distribution by a Boltzmann distribution, /| q ~ e~ E l T and take neutrinos to be either 
purely non-relativistic, i.e., E% = m s if p% < m s , or purely relativistic, i.e., E2 = P2 if P2 > m s . We define y = p/T, 
and integrate analytically Eq. (|2 . 2|) (with y constant) over temperatures T in the interval < T < Trh- We find // 
and fjj given by 



fi(T,y) 



s/y 



HT 



for < m s /y < Trh and 



e~ v dT 



(T,y) = J^l^-e-^dT 



s/y EiHl e ~^/T dT 

HT 
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(2.9) 



(2.10) 



for T < Trh < ms/y- Notice that due to the Boltzmann factor the contribution of the integrands at low temperatures 
is negligible. Thus, cutting the integrations at the decoupling temperature of active neutrinos or extending them to 
zero temperatures does not change the integrals in any significant way. 



2 We note that in this limit there is a factor of 2 with respect to Eq. fl'i.:>[ ) already present in the previous literature, which does not 
change our conclusions and we do not attempt to correct here. 

3 It was shown in Ref. Il0( that for masses lighter than 1 MeV this gives results which are correct within an order of magnitude [jj. 
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FIG. 1: Number density of (mostly) sterile neutrinos for 1 MeV < m s < 140 MeV and 5 MeV < Trh < m s . 



Once we have the distribution function, we then find the sterile neutrino number density as function of the ratio 
itIs/Trh and the temperature T, by integrating // and ///, 



r T( ms /T RH ) p 2 dp 

n s (m s I T RH , T) = / — T f H (p,T) 



o 



T(m s /T RH ) 



p 2 dp 

~2^ 



fi(p,T) 



(2.11) 



This procedure overestimates the final abundance, thus providing an upper bound on the actual abundance. In this 
way we find that the number density n s of heavy (mostly sterile) neutrinos for m s < 140 MeV is given by 



n s (x R H,T) ~ n a (T) sin 2 26 {^ff 2.1 x 1(T 3 e~ x 



7 6 + 144C(3)\ 2 3 x 7 f / 3 



3 4 C(3) 
7tt 4 
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24 + 144 C(3) + 12 x RH + - xln + 1 
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4tt 4 



(2.12) 



Here we have defined xrh 



"t-s/Trh- This number density is plotted in Fig. [TJ where n s /[n a sin 2 29 (m s /MeV) 3 ] is 
shown as function of xrh- Taking into account the subsequent decay of sterile neutrinos, the actual number density 
is n s (T) ~ n s (%Hi T) e^'/ 1 " 3 . In order to obtain analytical results for the cosmologically and astrophysically allowed 
regions in the parameter space (sin 2 28, m s ), in Section Hvl we will use the instant-decay approximation for t s smaller 
than the age of the Universe. 

We consider next the experimental bounds and then the cosmological and astrophysical bounds on the mass and 
mixing angle of heavy (mostly sterile) neutrinos. We will show that if the reheating temperature turns out to be 
sufficiently smaller than the neutrino mass, the cosmological bounds become irrelevant and the mixing angles for any 
given mass can be as large as the experimental bounds permit, making the detection at the reach of future laboratory 
experiments. 



III. EXPERIMENTAL BOUNDS 



In laboratory searches, no positive evidence of heavy (mostly sterile) neutrinos has been found so far in the mass 
range of interest, 1 MeV-140 MeV. Here, we review the most stringent bounds on the mixing angle with active 
neutrinos and show them in Figs. [2] and [3] (for further details see a comprehensive discussion in Ref. |22j|). 



5 



V - V mixin 

e s 




io" 8 io~ 7 io" 6 io" 5 io 4 io" 3 10" 2 10" 1 10° 

sin 2 29 

FIG. 2: Experimental bounds on the mass and mixing angle of sterile neutrinos mixed with i/ e in the mass range 1 MeV < 
m s < 140 MeV. The colored regions are excluded by searches of: i) kinks in /3-decays: horizontally hatched (red) area with 
short-dashed boundary; ii) sterile neutrino decays in visible particles: solid light gray (magenta) regions with dash-dotted 
contours; iii) peaks in the electron spectrum in pion and kaon decays: solid dark gray (blue) region delimited by a solid line; 
iv) neutrinoless double beta-decay: diagonally hatched (blue) area with long-dashed contour. See text for details. 



Let us consider first sterile neutrinos mixing with v e . For masses up to m s ~ 10 MeV, an important bound is provided 
by searches of kinks in the electron spectrum of /3-decays, which constrain the mixing angle to be sin 2 29 ^ 6 x 10~ 3 , 
as shown in the horizontally hatched (red) excluded region in Fig. O 

For higher masses, very robust bounds can be set by looking for additional peaks in the spectrum of electrons 
in leptonic decays of pions and kaons. The electron energy of the possible monochromatic line, given by Ei = 
{mir,K + ml — ml)/2m 7 c K , depends on the mass of the heavy sterile neutrino, while the branching ratio is proportional 
to sin 2 29. Here, m^ x is the mass of either the pion or the kaon respectively, m e is the mass of the electron and m B 
is the sterile neutrino mass. At present, bounds as stringent as sin 2 20 < 10~ 7 are obtained in this way (for a review 
see Ref. [23j]) which are shown in Fig. [2] as the solid dark gray (blue) excluded area. 

In neutrino-oscillation, fixed-target and collider experiments, if sterile neutrinos mix with active ones, a beam of v 8 
would be produced and would subsequently decay into visible particles. Assuming only charged current and neutral 
current interactions for v s , the absence of zArdecay signatures in past and present experiments allows one to put limits 
on the mixing term which controls the intensity of the v s beam and the decay time. A reanalysis of the Borexino 
Counting Test Facility and Bugey data yields sin 2 29 < 10~ 4 for m s < 10 MeV [24| at 90% confidence level (CL), while 
the data from the experiment PS191 [25] puts a bound which is strongly mass dependent, going from sin 2 29 < 4 x 10~ 4 
at m s ~ 12 MeV to sin 2 29 < 10~ 8 at m s ~ 300 MeV. These constraints exclude the solid light gray (light-magenta) 
regions in Fig. [51 

Finally, if sterile neutrinos are Majorana particles, they would contribute to the mediation of neutrinoless double- 
beta decay. The limit on the half- life time of this process can be translated into a bound on the mixing with v e , sin 2 9, 
which scales as m s for m s < 30 MeV and as m^ 1 for m s <; 400 MeV, with sin 2 29 < 5 x 10~ 8 at m s ~ 100 MeV. This 
limit excludes the diagonally hatched (blue) region shown in Fig. [5] 

Let us now consider the case of mixing with v^. As discussed above for the case of sterile neutrinos mixing with 



6 



v (v ) - v mixin 

\X v V s 




io" 8 io" 7 io" 6 io" 5 io" 4 io" 3 10" 2 IO" 1 10° 

sin 2 29 

FIG. 3: Experimental bounds on the mass and mixing angle of sterile neutrinos mixed with (or with v T ) in the mass 
range 1 MeV < m s < 140 MeV. For the case, the colored regions denote the bounds obtained from i) peak searches in the 
muon spectrum from pion and kaon decays which exclude the dark gray (blue) areas with solid contours; ii) searches of sterile 
neutrino decays excluding the light gray (light-magenta) region delimited by a dash-dotted line. The (cyan) long-dashed line 
represents the upper bound for sterile neutrinos mixed with v T . See text for further details. 



v e , peak searches provide very robust and stringent bounds on sterile neutrinos. Looking this time for peaks in the 
spectrum of muons in pion decays, it is possible to constrain the relevant mixing sin 2 29 to be typically < 10~ 4 in 
the mass range 4-34 MeV (for a detailed review see Fig. 1 of Ref. [III). Motivated by the KARMEN anomaly, which 
could be explained by the existence of a heavy (mostly sterile) neutrino with mass 33.9 MeV, very sensitive searches 
have been performed for this neutrino mass and the bound was found sin 2 29 < 9.2 x 10 -8 at 95% CL [27( (see the 
black horizontal line in Fig. [3]) . Similarly, sterile neutrinos with heavier masses can be probed in kaon decays and, in 
the mass range of interest, the derived bounds go from sin 29 < 4 x 1(T 4 at m s = 65 MeV to sin 29 < 25.2 (4) x 10 6 
for m s — 100 MeV (200 MeV) at 90% CL [28] . These peak-searches exclude the dark gray (blue) areas shown in 

By searching for the production of v s in pion and kaon decays and their subsequent decay, a reanalysis [261 ] of fixed- 
target data led to a stringent mass-dependent bound on sin 2 29, which exclude the light gray (light-magenta) area in 
Fig. [3] The bound is sin 2 29 < 10~ 4 for masses up to 35 MeV and reaches sin 2 29 < 3.4 x 10~ 7 at m s = 200 MeV, at 
90% CL. Similar bounds can be set on the product of the mixing angle with and with v e from the analysis of the 
same decay- searches data (a review is given in Ref. Hj|). 

Let us consider finally the case of heavy (mostly sterile) neutrinos mixed with v r . The only limits on these 
sterile neutrinos come from searches for v a decays. The most stringent bound is obtained by the reanalysis of data 
from the CHARM experiment, in which v s could be produced in D and r decays. However, in the mass range we are 
considering, this upper bound, shown in Fig. [3] by the (red) long-dashed line, is rather weak, such as sin 2 29 < 2 x 10~ 3 
for m s = 100 MeV. 
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IV. COSMOLOGICAL BOUNDS 



In what follows we obtain four different types of cosmological and astrophysical bounds on the parameter space of 
active-sterile neutrino mixing, (sin 2 20, m s ). Each of these limits is valid for a certain range of values of the heavy 
neutrino lifetime, t s . Whereas the first three bounds we describe have to do with the photons which are produced in 
the decay, the last one represents a bound on the sterile neutrino abundance at a particular epoch. Finally, we will 
also comment on bounds from core collapse supernovae observations. 

In the first place, we consider the diffuse extragalactic background radiation (DEBRA) spectrum and set bounds 
on the basis of not finding any unexpected result from heavy (mostly sterile) neutrino decay. Secondly, we obtain 
bounds using the non-observation of a distortion of the Cosmic Microwave Background (CMB) spectrum caused by the 
radiation from decay. Then, we calculate the limits based on the data from the primordial light clement abundances 
and how neutrino decay would affect them. Finally, we also use the upper bound from BBN on the extra number of 
relativistic degrees of freedom at that epoch to set a limit on active-sterile neutrino mixing. In our analysis, for t s < tjj 
with t\j ~ 14 Gyr the present age of the Universe, we use the instant-decay approximation, i.e., all decays happen at 
t = T a . This will allow us to obtain all the results analytically, while giving rise to accurate enough calculations for 
our purposes. 

For sterile neutrino masses in the range of interest, 1 MeV < m s < 140 MeV, the decay branching ratio of 1/2 
into e + e~ is about 10% (40%) for mixing with or v T (v e ). Due to inverse Compton scattering on the CMB 
with a Thompson cross section, the interaction length of these electrons is about 1 kpc (see for example Ref. or 
Fig. 5 of Ref. 30]) and the result of each of the interactions is a photon that shares a large portion of the incoming 
electron energy. Photons at these energies propagate for cosmological distances undisturbed by the CMB or infrared 
backgrounds and are subjected to different cosmological bounds, which depend on the time at which the photons were 
produced. 

Photons decouple from the plasma filling the Universe at the recombination epoch, t ICC ~ 1.3 x 10 13 sec. Thus, 
if the sterile neutrino decays happen after recombination, r s > t rec , the photons produced do not interact ever after 
and could leave an imprint in the DEBRA spectrum [31] . Such a signature has not been observed, and thus the 



photon flux must not be larger than the observed DEBRA energy flux 7 7 
the COMPTEL [32| instrument we obtain approximately 



EzdF 1 /dE 1 . From the observations by 



0.01- 



MeV 



cm- 1 sec sr 



for 1 MeV < / 7 to s < 30 MeV, where / 7 is the average fraction of the sterile neutrino mass that 
each decay. Thus / 7 m s is the average energy going into photons per decay. From EGRET data 



(4.1) 

oes into photons in 
we obtain 



Z 7 ~ 2.0 x 1(T 



MeV 



cnr sec sr 



(4.2) 



for 30 MeV < / 7 rn s < 140 MeV. In the instant-decay approximation, for t s < tu, the energy of the photons produced 
in the decays at t = t s redshifts from the time t = t s until today. For r s > tjj, because the decay rate increases with 
time for t < t s , most of the decays are happening at present, and we assume there is no significant redshift of the 
initial photon energy. For r s > t\j, we get 



LBm s 



47T 



< Iv 



(4.3) 



and for t rcc < t s < tjj, 



f~Bm s I — — 
1 tu 



2/3 



C T 
47T ' 



(4.4) 



where B is the branching ratio of the sterile neutrino decay into photons or charged particles. We label these bounds 
as "DEBRA" and it only appears in Fig. [7] 

On the other hand, the CMB radiation is emitted at recombination. Electromagnetic decay products produced 
sometime before recombination may distort the CMB spectrum |33l |34| . Non-thermal photons produced before the 
thermalization time ith — 10 6 sec are rapidly incorporated into the Planck spectrum. This happens through processes 
that change the number of photons, such as double Compton scattering ('ye — > 77e). If non-thermal photons are 
produced after i t h, i.e., if < t s < i rcc , the CMB Plank spectrum would be distorted. Current data pose very 
stringent upper bounds on possible distortions of this spectrum. For the earliest part of this last time interval, i.e., for 
tth — 10 6 sec < t s < 10 9 sec, photon number preserving processes, like elastic Compton scattering, are still efficient. 
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FIG. 4: Experimental and cosmological bounds as function of the mixing angle and mass of heavy (mostly sterile) neutrinos 
mixed with active ones for Trh = m s . The dark gray (blue) solid area and dark gray (blue) hatched area (for Majorana 
neutrinos) represent the experimentally excluded region for v E — v a mixing, the (red) short-dashed line and the (cyan) long- 
dashed line are the experimental upper bounds for — v s and v T — v a mixing, respectively. The cosmological and astrophysical 
bounds are indicated by the corresponding labels. Isolines for values of the heavy neutrino lifetime r equal to the start and end 
of the BBN epoch, the recombination time and the present age of the Universe are also shown. Only values of the reheating 
temperature larger than 4 MeV are considered. See text for details. 



These processes thermalize the photons not into a Plank spectrum but into a Bose- Einstein spectrum with a non-zero 
chemical potential fi. If the initial spectrum has fewer photons than a black body of the same total energy, then the 
chemical potential is positive \i > (if it has more photons, then fi < 0). So, for tth — 10 6 sec < r s < 10 9 sec, the 
energy released into photons in the decay for \fj,\ « 1 (the only values of [i allowed by observations) is 

^ ~ 0.714^ . (4.5) 

Pi 

The bound provided by the COBE satellite is < 0.9 x 10~ 4 at the 95% CL [35]. For later decays, 10 9 sec < t s < i rcc , 
the photon number preserving processes can no longer establish a Bose- Einstein spectrum. The energy released in this 
case is not thermalized but simply heats the electrons. Through inverse-Compton scattering this produces a distorted 
spectrum characterized by a parameter y, which for \y\ << 1 is related to the energy released in non-thermal photons 

as 



Ap 7 

Pi 



Ay . (4.6) 



The COBE upper bound on this parameter is \y\ < 1.5 x 10 5 [35[. In both cases, Eqs. (|4.5[) and (|4.6p . the fractional 
increase in the photon energy density due to the decay of the sterile neutrinos can be written as 

^El^^B^l (ZlV /2 < 6 x 10 -5 . (4.7) 
p 1 2.7 n-y Vsec/ 

All these bounds just described related to distortions in the CMB spectrum are labeled as "CMB" in Figs. [4j [5] and [7] 
For decays at earlier times, the best constraints on photons produced before the CMB thermalization epoch tth 
come from BBN which finishes by iendBBN — 10 4 sec. After £ e ndBBN, electromagnetic cascades can cause the photodis- 
sociation of D and 4 He. For photons produced in the time interval 10 4 sec < t s < 10 6 sec, the photodissociation of D 
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FIG. 5: Same as Fig. g] but for T RH = m a /3. 



poses the best limits. For earlier and later times, the overproduction of D due to the photodissociation of 4 He [33| far 
dominates its destruction, since the abundance of 4 He is about 10 4 times greater than that of D, and sets the most 
stringent constraints. The upper bounds taken from Fig. 3 of Ref. [H[ are 



VMeV/ 



n 7 \ r. 

for 10 4 sec < t s < 10 6 sec and m s > 2.2 MeV; 



5/2 

< io-a ( n^rr ) , (4.8) 



n s in _ 2 f 10 4 sec 



(4.9) 



for 10 6 sec < r s < 10 8 sec and m s > 28 MeV and 



VMeV/ \n 7 J V10 8 s 



1/4 

(4.10) 



for 10 8 sec < r s < 10 13 sec and m s > 28 MeV. We label these bounds as "Light Elements" in Figs. [H [5] and [7] 

For even earlier times, BBN data provides an upper bound on any source of extra energy density present in the 
Universe during the BBN epoch, f s tartBBN —0.1 sec < t < £ C ndBBN, as well as on extra radiation present during 
that period. The bounds are complicated in detail, but it is safe to say that if the extra energy density due to the 
presence of sterile neutrinos is very small, BBN will not be affected in any way. The bounds on extra contributions 
to the energy density during BBN are customarily presented in terms of the equivalent extra number of relativistic 
active neutrino species AN U . Thus, to be on the safe side, we simply require the AN V due to the presence of sterile 
neutrinos to be very small during BBN: 

AN V = ^~ <i. 4.11 

Pa 7T n a V MeV / V sec / 



We label this bound as "BBN" in Figs. HE] and 
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FIG. 6: Same as Fig. g] but for T RH = m s /W. 



Finally, there are also astro-physical disfavored regions in the (sin 2 29, m s ) space that we should mention. In principle, 
the energy loss into sterile neutrinos produced in core collapse supernovae explosions provides bounds on the mass 
and mixing angles of massive (mostly sterile) neutrinos. However, so much is not understood about the neutrino 
transport and flavor transformation in hot and dense nuclear matter, that conservatively the implicated m s ~ sin 2 29 
region can only be considered disfavored but not excluded. Only neutrinos with mass m s ^ 150 MeV could be 
emitted co piou sly in the supernova core and, for those masses, mixings sin 2 26 ) i; m = 3 x 10~ 10 is sin 2 29 is 10~ 2 are 
disfavored [IS, [36| . For small mixing angles the sterile neutrinos are not trapped within the supernova, thus they are 
emitted from the whole collapsing star, mostly from its core. The lower bound on the mixing angle is obtained by 
requiring that the sterile neutrino flux at Earth emitted by the supernova 1987A, F s , is not larger than the active 
flux, F a = 1.1 x 10 10 cm -2 , emitted by it. Thus, F s = F a for sin 2 29\\ m and for smaller angles 

sin 2 26» 

Fs= . 2o , F a . (4.12) 
sin 20iim 

This is relevant for the the last bound we will consider. The decay of heavy neutrinos emitted in supernova explosions 
would produce a flux of photons. The non-observation by the Solar Maximum Mission of any 7— ray counts in excess 
of the background for a time interval of i max = 223.2 sec after the arrival of the first 77 e 's from supernova 1987A, 
allows us to enlarge the disfavored regions of the parameter space. Following the analysis of Ref. [37| and neglecting 
the absorption of the photons produced in decays within the supernova, the 7 flux is given by 



dF 7 F s E, t / 7 B /■*— ^ e _ 2 Ei t/mt r + e -E-,IT (4 13) 



- 7 

dE 1 T V T 



Here T ~ 50 McV is the temperature at which sterile neutrinos are emitted from the supernova core. We use the 
3-er limits (AF 7 )3_ ff (£i,£j) on the 7 flux for the time interval considered obtained in Ref. [3711 for three energy 
bands (Ei,E f ), namely (4.1,4.4) MeV, (10,25) MeV and (25,100) MeV (see Table 1 of Ref. We calculate the 

disfavored region in the (sin 2 29, m s ) parameter space by requiring 



r E i dp 

J E ~d~E dEl < ( AF t)3-^, Ef) 



(4.14) 
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FIG. 7: Same as Fig. [4] but for a fixed value of the reheating temperature, Trh = 5 MeV and for heavy neutrino masses higher 
than 1 MeV. 



In the limit m s t >> 2 t max T, which holds in the region of interest, this condition can be written as 

G(Ei,E f ) < (AF 7 ) 3 _ CT ( J B i , E f ) , (4.15) 



with 



G{E i ,E f )^(e-^ T -e-^l T )+^ 



2 



(4.16) 

The band which gives the most restrictive limits for the range of sterile neutrino masses considered is (25, 100) MeV, 
although the band (10, 25) MeV gives comparable results. The two types of regions disfavored by core collapse 
supernovae arguments are the the dotted regions labeled "Supernova" in Figs. 01 [5l [5] and [JJ 

Figs. [4] to [7] show in the (sin 2 29, m s ) parameter space all the bounds on active-sterile neutrino mixing we have 
presented above, for T^m > 4 MeV and 1 MeV < m s < 140 MeV. For concreteness, in all the figures we have taken 
/ 7 = 1 and B — 0.1. Each cosmological bound is indicated by its corresponding label, as explained above. The regions 
globally excluded by laboratory measurements (Section IIII[) are: the dark gray (blue) area for v e — v s mixing, the 
region to the right of the red short-dashed line for — v s mixing and that to the right of the (magenta) long-dashed 
line for v T — v s mixing. In the case of v e — v s mixing the bounds for Majorana neutrinos are more restrictive than 
for Dirac neutrinos (see the hatched dark gray (blue) area labeled as "(f3/3)o v — decay"). Values of the heavy neutrino 
lifetime r equal to the relevant epochs in the history of the Universe are also shown. 

In Figs. |3J [5] and [SI we display the results for three different values of the ratio xb.h = to s /Trh> £rh = 1,3, 10, 
respectively. As expected, for Trh increasingly smaller than the sterile neutrino mass, the cosmological bounds become 
less restrictive, and when Trh < m s /10, the cosmological bounds become completely irrelevant and only experimental 
data are able to restrict the parameter space. This result can also be seen in Fig. [71 in which the reheating temperature 
is fixed to be T R h = 5 MeV. In this case all cosmological bounds become irrelevant for m s > 30 MeV. 
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V. CONCLUSIONS 

Sterile neutrinos are invoked in many extensions of the Standard Model of particle physics [HI, [38, 39, 40]. However, 
it is commonly assumed that the cosmological and astrophysical bounds on the mixings of sterile and active neutrinos 
restrict the range of their allowed values much more than laboratory data. In fact, sterile neutrinos with parameters 
suitable to be found in the near future in different experiments would have mixings with active neutrinos too large 
to be allowed by the standard cosmological assumptions about the pre-BBN era in the Universe, an era about which 
we do not have any observational information. The standard assumptions are few but very powerful: it is usually 
assumed that the temperature reached in the radiation dominated epoch before BBN was very high, that the Universe 
was radiation dominated then and that the entropy of radiation and matter is conserved. 

Here, we show that it is possible to evade most of the cosmological bounds by assuming that the temperature at the 
end of (the last episode of) inflation or entropy production, the so-called reheating temperature Trh, is low enough. 
We concentrate on massive (mostly sterile) neutrinos heavier than 1 MeV, having previously dealt with the lighter 
ones For low Trh, the production of sterile neutrinos is suppressed as shown in Fig. [1] For example, going from 
Trh = rn s to Trh = m s /10 there is a suppression of four orders of magnitude. We present the experimental bounds 
on sterile neutrino mixings and the cosmological bounds imposed by the diffuse extragalactic background radiation, 
the CMB, BBN and the abundance of light elements. We find that for Trh a few times smaller m s the cosmological 
bounds weaken significantly and for Trh = m s /10 they disappear completely. This is shown in Figs. [4j [5l [6j In 
Fig. [7] we keep instead the reheating temperature at a fixed value Trh = 5 MeV, close to the lower bound on Trh of 
4 MeV imposed by BBN and other cosmological data. In this case, no cosmological bounds remain for m a > 30 MeV, 
thus the only constraints on sterile neutrinos in this mass range come from terrestrial experiments. Hence, unlike 
in the standard cosmology, in low reheating temperature cosmologies it is possible to accommodate "visible" sterile 
neutrinos, i.e., sterile neutrinos which could soon be found in experiments. 

Cosmological scenarios with a very low reheating temperature are more complicated than the standard one. Al- 
though no consistent all-encompassing model of this nature exists at present, different aspects have been studied with 
interesting results, which suggest that a coherent scenario could be produced if an experimental indication would lead 
us to it. In fact, finding a particle, such as a "visible" sterile neutrino, whose existence would contradict the usual 
assumptions about the pre-BBN era, would give us not only invaluable information for particle physics, but also an 
indication of enormous relevance in cosmology: it would tell us that the usual assumptions must be modified, for 
example in the manner presented in this paper. 
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Appendix 



To be formally accurate, the derivation of Eq. (|2.2p should be carried out in the mass basis, since the distribution 
function is defined for species of definite energy. In the vacuum limit, in which we are interested, one can start with 
the Boltzmann equation in terms of the 2 by 2 density matrix p acting on the mass eigenstates vi and v> 2 , with 
elements p n , p 12 , p 2 \, P22 HH, 



if>= ~[H,p] -i{T,(p- p eq )} , 



(5.1) 



where p eq is given in terms of the equilibrium momentum distributions /° q for i = 1, 2, / 4 cq = (e^ Ei + 1) 1 {Ei 

and pi are the energy and chemical potential of the i mass eigenstate and without a lepton asymmetry pi = 0), 



Peq 



/r 

f 2 q 



(5.2) 



and the 2 by 2 matrix H is 



H = 



Ei 
E 2 



(5.3) 
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Neutrino production and destruction are represented by the anticommutator term in Eq. (|5.ip . which describes the 
coherence-breaking interactions. The matrix T can be written as 

I" cos 2 0(r/2 + £Ti) -cos6isin6l(r/2 + ^r) 1 

[ -cos0sin0(r/2 + £r) sin 2 6(T/2 + ST 2 ) J' { ' 

where T is the interaction rate for massless neutrinos, which is of order Gp (i.e., second order in the Fermi coupling 
constant), 5Ti are the corrections necessary when the incoming and outgoing particles have a non-zero mass TOj in 
the case of scattering, or when the two annihilating particles are of mass rrii and ST is the correction necessary when 
one of the two interacting particles has mass mi and the other has mass m,2- 

Replacing Eqs. (|5.2ll . (|5.3j) and (|5.4|) in Eq. (|5.ip we can find the equation for p'n, ^22, and p'12, which is the complex 
conjugate of p'21. Writing ip'12 = R + il and ip'21 = R — il, one can find the equations for the time derivative of the 
real and imaginary parts R and / in terms of R and I. As demonstrated in Ref. [42| . in the stationary point or static 
approximation, we can take R = I = which allows us to solve for R and /. Using the expression for R and / so 
obtained and assuming that the less massive eigenstate v\ maintains equilibrium, for small values of sin 2 28 one gets 
an expression for P22 which coincides with Eq. (|2.2[) in vacuum (once we take P22 = fz — fs and — /| q )- 
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